Mesoscopic superposition and sub-Planck scale structure in molecular wave packets 
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We demonstrate the possibility of realizing sub-Planck scale structures in the mesoscopic superpo- 
sition of molecular wave packets involving vibrational levels. The time evolution of the wave packet, 
taken here as the SU(2) coherent state of the Morse potential describing hydrogen iodide molecule, 
produces cat-like states, responsible for the above phenomenon. We investigate the phase space 
dynamics of the coherent state through the Wigner function approach and identify the interference 
phenomena behind the sub-Planck scale structures. The optimal parameter ranges are specified for 
observing these features. 
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Mesoscopic superposition of coherent states and their 
generalizations, such as cat-like states, have attracted 
considerable attention in the recent literature [j], 0> H|, 
since they show a host of non-classical behaviors. In a 
remarkable paper, Zurek Q demonstrated that appro- 
priate superposition of some of these states with a clas- 
sical action A can lead to sub-Planck scale structures in 
phase space. These sub-Planck scale structures in the 
phase space are characterized by an area h 2 /A. Apart 
from their counter intuitive nature and theoretical sig- 
nificance, the above scale has been shown to control the 
effectiveness of decoherence, a subject of tremendous cur- 
rent interest in the area of quantum computation and 
information. Zurek's realization made use of dynamical 
systems which exhibit chaotic behavior in the classical 
domain. Recently a cavity QED realization involving the 
mesoscopic superposition of the compass states have been 
given Q. In principle, one could also use superpositions 
of cat-like states arising in quantum optical systems with 
large Kerr non-linearity ■ 

In this paper, we demonstrate the possibility of real- 
izing sub-Planck scale structures in the mesoscopic su- 
perposition of molecular wave packets, which involves vi- 
brational levels. The time evolution of an initial wave 
packet, taken here as the SU(2) coherent state (CS) of 
the Morse potential produces cat-like states. These arise 
due to the quadratic dependence of the energy on the vi- 
brational quantum number. The superposition of these 
states is responsible for the above phenomena. We study 
the spatio-temporal structure of these states, paying spe- 
cial attention to the fractional revival, which gives rise 
to four coherent states required for the observation of 
the sub-Planck structure. This structure can be clearly 
explained through the interference phenomena in phase 
space. For this, we investigate the phase space dynamics 
of the coherent state through the Wigner function ap- 
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proach and identify the optimal parameter ranges for a 
clear observation of these features. 

Morse potential is well-known to capture the vibra- 
tional dynamics of a number of diatomic molecules 
[1, 0, H, H E3- It is worth mentioning that the phe- 
nomena of revival and fractional revival Jll|, [3 have 
been experimentally observed in wave packets involving 
vibrational levels [lj] . Creation of the wave packets and 
observation of their dynamics are carried out through 



pump-probe method [15j. The control and analysis of 



molecular dynamics is achieved through ultrashort femto- 
second laser pulses [l6| ■ Fractional revival can be probed 
by random-phase fluorescence interferometry [17J. Re- 
cently, cat-like states, arising in the temporal evolution 
of the Morse system, have been proposed for use in the 
quantum logic operations [Tc| . 

Morse potential describing the vibrational motion of a 
diatomic molecule has the form 



V(x) = D(e 



-28x 



2e~ f3x ) 



(1) 



where x = r/ro — 1, rn is the equilibrium value of 
the inter-nuclear distance r, D is the dissociation en- 
ergy and (3 is a range parameter. We will be consid- 
ering HI molecule, as an example, which has 30 bound 
states, with (3 — 2.07932, reduced mass p — 1819.99 a.u., 
r = 3.04159 a.u. and D = 0.1125 a.u. Defining 
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and 



p 2 h 2 



E, 



(2) 



eigen functions of the Morse potential can be written as 

^(0 = Ne-^ s/2 K(0, (3) 

where £ = 2\e~ t3x \ < £ < oo, and n = 0, 1, [A - 
1/2], with [p] denoting the largest integer smaller than 
p, so that the total number of bound states is [A — 1/2]. 
The parameters A and s satisfy the constraint condition 
s + 2n = 2A - 1. 

Note that A is potential dependent, s is related to en- 
ergy E and, by definition, A > 0, s > 0. In Eq. ([3]), 
L s n (y) is the associated Laguerre polynomial and N is 
the normalization constant: 
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FIG. 1: (Color online) \d m \ 2 plotted as a function of m for 
the Morse potential of HI molecule for different values of a. 
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Quite some time back, Nieto and Simmons gave a 
minimum uncertainty coherent state for Morse oscilla- 
tor considering suitable conjugate variables fioj ] . Later, 
Benedict and Molnar [2(| also found the same CS through 
super symmetric quantum mechanical method. This was 
used to describe the cat states of the NO molecule (2lj . 
This CS involves infinite number of bound states, not be- 
longing to the same potential [22| . Morse potential has 
a finite number of bound states. Hence it is natural to 
expect an underlying SU(2) algebra. Recently, Dong et 
al, [23[ have obtained the SU(2) generators J+, J_ and 
J which satisfy the algebra at the level of wave function 



J+,J_ 



V£(0 = 2J ^(C). 



(5) 



The SU(2) Perelomov coherent state of the Morse sys- 
tem is obtained by operating the displacement opera- 
tor e aJ +~ a J - on the highest bound state n', defined by 
J+ipn'iO ~ 0- Using disentanglement theorem, the co- 
herent state modulo normalization becomes 



x(0 



,-aJ-,/,X 



^, - ay/n'(s + n' + l) + 

(-")" ' 
(n')\ 



+ \' y/n'\{s + n' + l)(s + n' + 2) 



xy/....(s + 2n') % 



(6) 



As is explicitly seen the above CS involves only the 
bound states, which are finite in number. This is due 
to the fact that the underlying group here is a compact 
group [24) | . For the purpose of our analysis, we consider 



this wave packet. We have checked that, superposition of 
Morse eigen states with suitable Gaussian weight factors, 
also reproduces the sub-Planck scale structure. 

Simplifying the above expression, we can write it in a 
compact form: 



where 



dm — 



x(0 



(- a y 



dm 



m=0 



(n' — m)! 
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Fig. [T] shows |e? m | 2 distribution of HI molecule for var- 
ious values of a. For smaller values of a, |d m | 2 is peaked 
at higher values of m, where the anharmonicity is larger. 
The corresponding initial CS wave packet is not well lo- 
calized and has an oscillatory tail. With the increase of a, 
\d m \ 2 distribution moves towards the lower levels and the 
wave packet's oscillatory tail gradually disappears. For 
larger values of a, only the lower levels contribute to form 
the CS wave packet, where the effect of anharmonicity is 
rather small. Hence, it is clear that the choice of the dis- 
tribution is quite crucial in the wave packet localization 
and its subsequent dynamics. 

Temporal evolution of CS state wave packet is given 

by 



(9) 



with E m = -(D/A 2 )(A 



1/2) 2 . This quadratic en- 



ergy spectrum yields classical and the revival times given 
by T d = T rcv /(2A - 1) and T rcv = 2n\ 2 /D respectively. 
More interestingly, when t takes the values 7 -T ICV , where 
r and q are mutually prime integers, the CS wave packet 
can be written as a sum of classical CS wave packets : 

l-i 

x(U) = J2 a p ( r /« t ^-pI 1 T d)]> ( 10 ) 



where 

Xd(e,*) =53 d m#UOe*p[-27rimi/T cl ]. (11) 

m 

Amplitudes are determined by 
1 

a p = — exp \2ni(m 2 r/q — mp/Z)] , (12) 
m 

where Z = q/2 when q is an integer multiple of 4 and 
Z — q, in all other cases. 

Fig.[5]shows the CS wave packet in the co-ordinate rep- 
resentation, where the revival behaviors at T rev /4 and 
T rev /8 are not transparent. We will now clarify the 
phase space picture of the wave packet at fractional re- 
vival times by using the Wigner function approach. We 
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FIG. 2: (Color online) Wave packet distribution in coordinate space for HI molecule, where a — 1.4, /3 = 2.07932. Plotted here 
is lx(£,*)| 2 as a function of x (where £ = 2A exp[-/3x]) for (a) t=0, (b) t = T rcv /8, (c) t = T rcv /4 and (d) t = T rcv /2. 



will also show that the interference phenomenon in phase 
space involving the cat-like states gives rise to the sub- 
Planck scale structure. 

The Wigner function can be written as 



W(x,p,t) = / x*(x-x',t) 

xx(x + x',t)e- 2lpx ' /h dx' , (13) 



where x is the scaled co-ordinate and p is the correspond- 
ing scaled momentum and also x(x) = x(£)- 

Wigner functions at instances of fractional revival can 
be explained by making use of the decomposition of 
Eq. (fTU]) . At t = Xicv/8, for example, the CS wave packet 
splits into four classical wave packets: 



/ j. 2rev -*el \ 

Xd(t, -g -£) 



Xc/(f,-5 tt) 



+ XcM,—^ ^-)J- 



(14) 



Defining 

(even. odd) 



m ovon „dd 



dmipmiO exp[-27rim— ] 

J-r] 



(15) 

expression Eq. (jT4j) can be rewritten in a simpler form: 



(16) 

The above expression plays a crucial role in the ex- 
planation of the phase space behavior at T rev /8. Sub- 
stituting this in Eq. (fTB")) . the Wigner function at t = 
T rev /8 can be written down as a sum of three terms: 
W(x,p,T Iev /8) = W^ cvcn ^ + + W( int \ where 

Ty(ovcn) an( j T^/(odd) are j_ ne Signer functions correspond- 
ing to the first and second terms on the right hand side of 



Eq. (fl"6f and T4 y ^ nt - ) is the contribution from the interfer- 
ence between these two terms. In Fig. [31 we have plotted 
W(x,p,T Tev /8) and its constituent parts for two values 
of a. Note that both W^ ovcn ) and W^( odd ) are Wigner 
functions of cat-like states. Each consists of two distinct 
peaks corresponding to two mesoscopic wave packets, 
and an oscillatory structure at the middle due to quan- 
tum interference between them. Furthermore, VF( cvcn ) is 
along the east-west direction whereas W^ odd ^ is along the 
north-south. This is because the time arguments of Xc™ 1 
and x°f d differ by T c i/4 in Eq. fT5)). The superposition 
of the interference regions of VF~( cvcn ) and W^ od<i ^ gives 
rise to the sub-Planck structure in Fig. [31(d) . It is worth 
pointing out that W lnt , as plotted in Fig. [3(c), gives the 
off diagonal interferences of compass-like states produced 
at T rcv /8. 

As seen in Fig. [H for higher values of a, the initial 
wave packet involves lower vibrational levels for which 
the turning points are nearer, resulting in a decrease in 
the span of the phase space variables. In this case, the 
area of overlap between the two interference structure 
increases and the number of ripples become less. As a 
consequence, the sub-Planck scale structure at the middle 
becomes more prominent as shown in the bottom array of 
Fig. [3j The four mini-wave packets, produced at T rev /8, 
are not equi-spaced and not of same size. The asymmet- 
rical nature of the Morse potential is the main reason 
behind this. We also analyze numerically the expecta- 
tion values of position and momentum at t — T rov /8 for 
different values of a. The uncertainty product (AxAp), 
obtained from this analysis, is 5.5914 for a — 1.4 and 
2.56404 for a — 2.5 in the unit of h = 1. The classical 
action is defined by A w AxAp and the corresponding 
dimension of the sub-Planck scale structure is a « h 2 /A 
which easily comes out to be 0.179 for a — 1.4 and 
0.39 for a — 2.5 respectively, implying the sub-Planck 
scale structure. Note that for smaller values of a the 
area becomes more sub-Planck. 

In conclusion, we demonstrate that the interesting sub- 
Planck structure in mesoscopic quantum systems can in- 
deed be realized in the temporal evolution of vibrational 
wave packets. This is clearly present, where four wave 
packets are produced in the temporal evolution. The co- 
herence parameter a plays a crucial role in the formation 
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FIG. 3: (Color online) The Wigner function W(x,p, t) and its constituent parts at t = T rov /8 as a function of x and p for a = 1.4 
(top row) and a = 2.5 (bottom row). Shown here are the contour plots of (a)W (cvon) ; (b)W (odd) ; (c)W (int) and (d)W(x,p,t). 



of this structure; one needs the low-lying states for a clear ence through coupling to rotational and other vibrational 
observation of this structure. It is worth pointing out levels. Recently methods like closed-loop control 25j has 
that, the vibrational wave packets are prone to decoher- been devised to minimize the decoherence effect. 
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